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Abstract 

Following arXiv: 1012.2107 we show that in global de Sitter space its isometry is broken by 
the loop IR divergences for any invariant vacuum state of the massive scalars. We derive 
kinetic equation in global de Sitter space. It follows from the Dyson- Schwinger equation of the 
Schwinger-Keldysh diagrammatic technique in IR limit and allows to understand the physical 
meaning and consequences of the loop IR divergences. In many respects the isometry breaking 
in global dS is similar to the one in the contracting Poincare patch of de Sitter space. Hence, 
as a warm up exercise we study the kinetic equation and properties of its solutions in the 
expanding and contracting Poincare patches of de Sitter space. Quite unexpectedly we find 
that under some initial conditions there is an explosive production of massive particles in the 
expanding Poincare patch. 



I. INTRODUCTION 



De Sitter's (dS) space isometry is crucial for its stability. Obviously if it is not broken, 
then the cosmological constant can not be secularly screened. There are states of the massive 
scalar field theories in which they respect dS isometry at tree level . But in the loops the 
situation becomes more complicated. One loop calculation in the expanding Poincare patch 
(EPP) of dS space shows that the isometry is broken for all a-vacua except Bunch-Davies (BD) 
one (see also the Appendix of {4]). 
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However, we find it rather unphysical to address the question of the stability of the system 
under such conditions when all its symmetries are respected exactly We propose to consider 
excitations above the highly symmetric state and to trace where do they evolve in the future 
infinity. For the finite particle density on top of any a-vacuum (including BD one) dS isometry 
is broken even at tree-level The obvious question to ask at this point is whether this 

density decays or rapidly increases with time? Rephrasing this, we would like to understand if 
the system has the dS invariant ground state. 

The standard answer to this question is that the density will decay due to the rapid ex- 
pansion of the space. But the situation happens to be more complicated and actually quite 
counterintuitive. We show in this paper that the proper quantity to study is the number density 
per comoving volume. (Unlike the physical volume, which inflates with time, comoving volume 
remains constant.) Note that the number density per comoving volume of the spatially homo- 
geneous dust remains constant throughout the history of dS space independently of whether 
the space is expanding or contracting. But such a density will grow unboundedly if one will 
turn on only the process of particle production by the gravitational field. That is, in particular, 
the source of the loop IR divergences and large IR contributions, which are under study in this 
paper. Of cause for the issue of the backreaction it is important to know how fast is such a 
growth in comparison with the expansion of the space, especially when the loss processes are 
also taken into account (apart from the production ones)*. The main goal of this paper is to 
address this problem. 

In global dS the situation is more intricate. There is one-loop IR divergence even for massive 
scalars [5|. It reveals itself via impossibility to take the moment of turning on selfinteractions, 
t , to past infinity. We are going to see that t plays the role of IR cutoff. Thus, unlike the 
situation with the BD initial state in EPP, in global dS isometry is broken in the loops for 
all possible vacua (including Euclidian one). Furthermore, the presence of such IR divergences 
means that the grows of particle density wins against the expansion of the space-time. 

However, it remains unclear whether after the summation of all loop leading IR contributions 
one can take to — > — oo or cannot. We show that this problem is related to the one of the 
presence of a stationary state in dS QFT. If such a state does exist, then t can be taken to 



* Note that turning on only single process of the particle kinetics, rather than all of them, may cause an 
instability (see [3] for the detailed discussion). However, as we show in many cases these instabilities do 
disappear when one takes into account all processes. But then new ones of different kinds do appear. 
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past infinity and the dS invariance is restored. Otherwise the dS isometry is not restored at 
higher loops and the backreaction due to the particle production cannot be neglected. It was 
explained in 0, Q that to approach such a problem one has to solve the kinetic equation, which 
is the tool for the summation of the leading IR contributions in all loops. Below we derive this 
equation in global dS. 

Although in EPP there are no IR divergences, one has there large IR contributions. As the 
result, one can derive the kinetic equation in EPP (4, 6]. In the latter papers the problem of 
the existence of the stationary state in EPP was addressed. The conclusion was that it can 
exist with the appropriate initial conditions. 

In the contracting Poincare patch (CPP) the situation happens to be similar to the global dS. 
There, unlike EPP, we also encounter loop IR divergences. In the main body of this paper we 
explain the reason of such a difference between EPP and CPP. We derive the kinetic equation 
in CPP. Such equations in EPP and CPP can be simply mapped to each other, because we are 
dealing with the spatially homogeneous states and study the density per comoving volume. As 
the result there is a quite simple map between the solutions of the kinetic equations in EPP 
and CPP. In particular, we find that there is the stationary solution even in the CPP, which is 
related to the above mentioned stationary solution in EPP. Its presence means that under the 
appropriate initial conditions and in the circumstances of the perfect spatial homogeneity, dS 
isometry can be restored, after the summation of all leading loop contributions, even in CPP. 

But one of the main messages of this paper is that there is also a solution of the kinetic 
equation (both in CPP and EPP) which describes explosion, within a finite proper time, of 
the particle density. The presence of such a solution means that under some initial conditions 
backreaction from the particle production can become enormous even in EPP. The last type of 
the situation should lead to the deformation of the dS metric and probably bring the system 
to the FRW space-time with secularly decreasing cosmological constant. 



A. Subtleties 



To proceed let us point out some subtleties with field theory in curved space-times, in general, 
and in dS space, in particular. For illustrative reasons consider 4D massless conformaly coupled 
scalar field with the quartic selfinteraction, [□ + R/6] = A0 3 , on global dS background. We 
are interested in the following classical scattering problem. Selfinteraction is adiabatically 
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turned on somewere in the far past and then eventually switched off in the far future. At 
past infinity we excite a single free wave and would like to ask the following question: will the 
solution of the equation under consideration (at the linear order in A) contain three free waves 
in future infinity or will not? 

In the metric with noncompact spatial slices ds 2 = Jj- [drj 2 — dx 2 ], 77 G (—00, +00) (see fig. 
1), the problem can be mapped, via conformal transformation, to the one in Minkowski space. 
The intensity of the process in question, at the leading A order, will be proportional to the 
5-function ensuring energy conservation |2J. The latter will put an obstruction to the process 
under consideration. 

If, however, one will address the same problem in the metric ds 2 = dt 2 — cosh 2 (t)dVl 2 = 
cos 2( r ) [dT 2 — dfl 2 } (see fig. 2), then it can be converted to the one in Einstein static universe 
with the compact range of time, T 6 [— | , +§]. As a result instead of the aforementioned 8- 
function one will encounter sm ^^ , where AE is the energy deficit in the process in question. 
Hence, the process will be allowed. 

The resolution of the apparent disagreement with the general covariance is that the result of 
this scattering process strongly depends on the initial/boundary conditions. Approaching the 
same problem with the use of different ways of spatial slicing of dS space we implicitly assumed 
different initial conditions. In fact, while in the case of noncompact spatial sections one has 
such harmonics as e~ iWV , for the situation with the compact sections the harmonics are e~ ttT . 
The map between these harmonics under the coordinate and conformal transformations should 
be of the main concern for those who are studying the scattering problem under consideration. 

The observations of the above paragraphs are relevant for us because below we are going 

to study IR contributions, which of cause are sensitive to the conditions at the boundaries of 

various patches of the entire space-time. Also they give one a hint for the explanation of the 

disagreement between the conclusions of [8fl and those of the present paper. The spatial slicing 

fl 

of fig. 1 was used in [8j|. Although in the latter paper the entire dS space was considered, 
the BD state was imposed at the boundary separating contracting CPP and EPP. We do not 
quite understand the physical meaning of this approach, but we can see that in such a way it 
is possible to respect dS isometry at each step of the construction of the correlation functions. 
In fact, as we have already pointed out above, to respect dS isometry in all loops one has to 
impose BD state exactly at the boundary separating EPP from CPP. 

At the same time, in the kinetic theory approach one can use both spatial slicings of fig. 1 



Figure 1: Penrose diagram of dS space. Left and right sides of the depicted rectangle are glued to 
eachother. Here we show Cauchy surfaces (e.g. constant rj slices) for the global dS metric of the form 
ds 2 = -y [dr] 2 — dx 2 ] , where r\ G (—00, +00). 



Figure 2: Here we show Cauchy surfaces (e.g. constant t or T slices) for the global dS metric of the 
form ds 2 = dt 2 - cosh 2 (t)dn 2 = [dT 2 - dfl 2 ] , where t G (-00, +00) and T G [-§, +§ ]. 

and fig. 2. In this setting there is nothing behind the initial Cauchy surface and the dS isometry 
is broken, at least in the loops, unless we start with the exact BD state at the boundary of 
EPP. 
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B. The setup of the problem 

Our eventual goal is to see whether the back-reaction from the particle production in dS 
space is negligible or not. But in this paper we neglect back-reaction and consider self- 
interacting massive QFT on the fixed dS background. We restrict our attention to the behavior 
of the solutions of the kinetic equation in the IR limit and try to see if the assumption of the 
negligible backreaction is self consistent. Note that in EPP (CPP) the IR limit of the physical 
momentum corresponds to future (past) infinity. 

We study D- dimensional minimally coupled real scalar field theory: 



Dots here stand for the higher power self-interaction terms, which make the theory stable. 
However, below we are going to consider formulas only due to the unstable cubic part of the 
potential. The reason for that is just to simplify equations. This instability does not affect our 

n 

conclusions, because similar IR contributions do also appear in massive 4 theory [9j] or in the 
other nonconformal theories in dS space 10], 12]. jl^j], jl^ ]. 

To set the notations we describe here some of the properties of this space and of free scalar 
harmonics on it. It is the hyperboloid Xq — X 2 = —1, i — 1, . . . , D in the (D + l)-dimensional 
Minkowski space ds 2 = cLXq — dXf. Throughout this paper we set the curvature of the hyper- 
boloid to one. The dS isometry group is SO(D, 1) Lorentz group of the ambient Minkowski 
space in question. 

The global dS metric is induced through the solution of the equation Xq — Xf = — 1 via 
Xq = sinht and Xj = Ui cosht, where u 2 = 1 defines the coordinates on the (D— l)-dimensional 
sphere. The induced metric is then ds 2 = dt 2 — cosh 2 tdfl 2 . 

The EPP is the halfspace of this hyperboloid. It is defined by the condition Xo > Xp. 
Any other EPP is obtained by a rotation of the latter patch around X axis, i.e. there is a 
family of EEPs, which is parametrized by the angle of the rotation around Xo axis. Such a 
patch of dS respects only some subgroup of the whole isometry group. Its induced metric is 
ds 2 = dt 2 — e 2t dx 2 = -5 [drf — dx 2 }, where 77 = e~* is the conformal time, which is flowing in 
the reverse direction — from infinity (in the past) to zero (in the future). 

The CPP also is a halfspace of the hyperboloid, which is defined e.g. by the condition 
Xq < Xd, i.e. it is complementary to EPP within the global dS. The induced metric on the 



7 



CPP is ds 2 = dt 2 — e~ 2t dx 2 = [d-rf — dx 2 ], where r] = e t G [0, +00). 

For large enough D the theory ([1]) becomes non-renormalizable. However, as we explain in 
the concluding section, if there is a stationary state in this theory, then in the IR limit we do 
not care about UV divergences and renormalizability of the theory in question. We assume that 
all couplings in all equations take their physical values, i.e. all UV divergent (~ A 2 log A) or 
finite (~ A 2 ) contributions are absorbed into their renormalization. It is straightforward to do 
that during the derivation of the kinetic equation from the Dyson-Schwinger one. At this point 
one should keep in mind that for the propagators which have proper Hadamard behavior the 
UV divergences in dS space are the same as in flat one. Only BD harmonics provide the proper 
Hadamard behavior in EPP and CPP. In global dS the same is true only for the Euclidian 
harmonics. 

The scalar field harmonics in global dS satisfy the Klein-Gordon equation: 



% + {D-2) tanh(t) d t + m 2 - 



(2) 



coslr(t) 

where Afl_i(fi) is the Laplacian on the (D — l)-dimensional sphere. Any solution of this 
equation can be expanded in harmonics as follows <f) p (t,Q) = g p (t) Y pr ^(Q), where 5^,^,(0) are 
(D — 1) -dimensional spherical functions and g p (t) can be expressed through the hypergeometric 
function l|> 3- 

For large enough p we can somewhat simplify eq.(j2]) in the limit t — > ±00. In fact, under these 
conditions and after the change of variables T]± = e Tt the global dS metric can be approximated 
by ds\ ~ [df] 2 - — dx 2 ~\ . Here conformal time 77+ flows towards zero (— > +0) in future infinity 
and, thus, covers a part of the EPP; similarly rj_ starts from zero (OH >) at past infinity, and 

D-l 

covers a part of the CPP. Furthermore, the substitution of <p(t — > ±00) w rj± 2 h±(prj±) e Ttpx 
into the simplified (by the limit t — > ±00) form of (J2J) reduces it to that of Bessel for h±(prf). 
The index of the obtained Bessel equation is ifi, where \x = J m 2 — (^j^~) 2 - 

Throughout this paper we are going to work with the massive fields from the principal series, 
m > (D — l)/2. It can be shown [9(, however, that for the fields from the complementary series, 
m < (D— 1)/2, one- loop IR divergences have the similar logarithmic character to that of massive 
fields. However, for the fields from the complementary series it is not clear how to write the 
kinetic equation of the conventional form, because for them the necessary separation of scales 
between the time dependence of the harmonics and of the particle density is not obvious. 
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II. CONTRACTING POINCARE PATCH 



The reason why we start our discussion with CPP is that the physics of IR divergences there 
is relevant for the field dynamics in global dS. 

To study loop dynamics in the non-stationary situation one has to use the in-in or Schwinger- 
Keldysh diagrammatic technique. Feynman rules in EPP can be found in [if]]. In CPP the 
rules are very similar. 

Each particle in the non-stationary situation is described by the matrix of propagators. 
The constituents of this matrix are the retarded and advanced propagators D R ' A , which carry 
information about the particle spectrum of the theory in question, and the Keldysh propagator 
D K , which carries the statistical information, i.e. shows which levels are populated. 

Due to spatial homogeneity of CPP together with the initial state, we find it con- 
venient to perform the Fourier transform of all quantities along the spatial direction 
x: Dp' R ' A (r]i,ri2) = / d D ~ 1 x e ipx D K > R > (7)1, x; 772, 0). Then at tree level the Keldysh 
propagator is 1)^,(771,772) = (^i^) -2 ^ Re [h(pT)i)h*(pr)2)], retarded and advanced prop- 
agators are £^(771,772) = 6 (r]i - V2) 2 (r)^)^ 1 Im [h(prj 1 )h*(prj 2 )} and 1)^(771,772) = 
—d(rj 2 — rii) 2 (r/i^)^ 1 lm[h(pr)i)h*(prj2)}; h{prf) is a solution of the Bessel equation, which 
we do not specify so far. The choice of the solution defines which harmonic basis we are using 
to quantize the theory. In other words, it specifies what we have chosen as the background 
state in our problem. 

In the limit 77 = s/rji^ — > 00, 771/772 = const and 770 — > the leading one-loop IR contribution 
to D*f (7)1,7)2) is contained within the following expression 



D i P (Vh m) = (Vi V2) 2 MPVi) ^*(P^2) n p (r)) + h(pr)x) h{pq 2 ) k p (v) + c - c -\ 



where n p (rj) 



A 2 Sd-2 



d D l q I I dr} Z dr\± (773774) 



and ftp (77) 



(2vr) D - 1 

Xh* (pr] 3 ) h (pr] 4 ) h*(qr] 3 ) h(qn A ) h*(\q - ^773) h(\p - q\r] 4 ), 
2\ 2 S D ^ 



d q / dr/ 3 / ^774(773774 



D-i 
2 X 



//n 



(2tt) d - 1 

xh* (pr) 3 ) h* (pn A ) h* (q7] 3 ) h(qn 4 ) h* (\q - p\n 3 ) h(\q - p\r] A ) . (3) 

The details of a similar calculation in EPP can be found in (4-6]. The largest IR contribution 
to (J3j) comes from the region of integration over dq, where q ^> p. It is instructive to study the 
character of such divergences for different choices of the harmonics h(pn). If h(x) oc T-L^(x) 
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the IR divergence is present both in n p and n p and is proportional to log (77/770) = t — t , if 
prj <C 1, which is the proper time elapsed from 770 = e*° to 77 = e*. If, however, prj ^> 1, then 
the divergence is proportional to log (l/p77o). The physical motivation for the latter behavior 
is given at the end of the next subsection. 

The in-harmonics of the CPP are proportional to the Bessel function Y ifJL (x). In this case 
77^(77) diverges as log (77 / 770) , if pi] ^ I5 an d as log (l/p?7o), when prj 3> 1. At the same time, 
ftp (77) does not have any IR divergence. We interpret this fact as that the backreaction, due to 
the particle production on the background state, corresponding to the in-harmonics, is small. 
That is the reason why we choose these harmonics for the proper definition of the quasi-particles 

in the IR limit |6| (see also the discussion below). For any other type of a-harmonics the IR 

(2) 

behavior of n p and K p is similar to that of Ti\ (pf]) harmonics. 



A. The physical interpretation of the one-loop IR divergence 

Now we are going to show that n p given in ([3]) is the total particle number density per 
comoving volume, {a p a p ) , created by the gravitational field during the time passed from 770 to 
77. Similarly k p is the measure of the strength of the backreaction, (a p a_ p ), due to the particle 
production process on the background state during the same time period. 

To explain the statement of the last paragraph, let us forget for the moment about the 
presence of k p and consider the situation in flat space with the plane wave harmonics. If the 
time of the turning on selfinteractions is to and the time of the observation is t, then the kinetic 
equation is (the derivation can be found e.g. in {(]]): 



dn r 



dt' cos [(-e(p) + e(q) + e(p - q)) (t - t')] (1 + n p ) n q 77 |p _ a| - n p (1 + n q ) (1 + n [p ^ 



dt 

'\p-i\ 



oc -A z 



d D ~\ 



:(p) e(q) e(p - q) 



(t) 



+2 / dt' cos [(-e(q) + e(q - p) + e(p)) (t - t')\ n q (1 + n\ q - p \) (1 + n p ) - (1 + n q ) n\ q ^ p \ n p (t) 



to 



+ / dt' cos [(e(q)+e(q+p)+e(p)) {t-t')\ (1 + n q ) (1 + n q+p ) (1 + n p ) -n q n q+p n p (t) W4) 



to 



where e(p) is the energy of the particle with the momentum p. The RHS of this equation is the 
collision integral Icon- I n the limit (t — to) — > 00 the dt' integrals transform into the (^-functions 
ensuring energy conservation in some processes, which we are going to define now. 
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The expressions multiplying the three dt' integrals in (j3J) have the following meaning. The 
first term describes the competition between the following two processes. One of them is that 
the particle p decays into two excitations — q and p — q. It corresponds to the term n p (1 + 
n q ) (1 + n\p- q \) and appears with the minus sign in the collision integral, because it describes 
the loss of the excitation with the momentum p. The inverse gain process, corresponding to 
the term (1 + n p ) n q n\ p - q \ with the plus sign, is such that two particles, q and p — q, are merged 
together to create the excitation p. 

The second term in (j3J) also describes two competing processes. The first of them is such 
that the particle p joins together with another one, q — p, to create the excitation with the 
momentum q. This is the loss process. The inverse gain process is the one in which the particle 
q decays into two, one of which is with the momentum p. The coefficient 2 in front of this term 
is just the combinatoric factor. 

Similarly the third term describes two processes. The gain process is when three particles, 
one of which is with the momentum p, are created by an external field, if any. The loss process 
is when three such excitations are annihilated into the vacuum. 

All these six types of processes are not allowed by the energy-momentum conservation for 
massive fields with <p 3 selfinteraction in flat space-time. Hence, the collision integral is vanishing 
as {t — t ) °°- 



However, in dS space they are allowed 



17|]-20|, because there is no energy conservation. 



As we show below the collision integral in dS space can be obtained from (j4]) as follows. One 
exchanges dt 1 for drj', multiplies by the proper weight coming from the metric, and uses under 
the corresponding integrals the dS harmonics, g p (rj) = rf~^ h(pr}), instead of the plane waves, 
e -ie(p)t J e^pj. Then, in place of the above mentioned 5-functions, which ensure energy conser- 



vation, there are some expressions whose physical meaning is that they define the differential 
rates (per given range of q) of the processes described in the above three paragraphs. 

Now put n p to zero under the collision integral (jlj) with the substitutions described in 
the previous paragraph. The only term which survives is present in the last line of (jl]) and 
describes particle creation by the gravitational field. If one takes the integral of this term over 
the conformal time from r) to rj he obtains exactly the expression for n p (rj) represented in fl3]). 
Using the collision integral for k p , derived in {4, ^ (see also the discussion below), one can make 
similar observation about the origin of the corresponding expression in ([3]). 

Thus, in the non-stationary situations, when Icon is not zero, it is natural to expect the linear 
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divergence, oc (t— 1 ), in n p and K p . That is, indeed, should be true for the homogeneous in time 
background fields (such as e.g. constant electric fields in QED), when particle production rates 
are constant in time. But, in time dependent metric of dS space one encounters a bit different 
situation. In particular, in EPP the largest IR contribution has the form of A 2 log{prj) and there 
is no IR divergence as one takes the moment of turning on selfinteractions to past infinity of the 
EPP (?7o — > +00). That happens because the creation of particles with comoving momentum 
p effectively starts at 77* ~ l/p rather than at r]o — >■ +00. Similarly in the case of CPP the 
particle creation process goes on until 77* ~ l/p and then stops. We encounter such a situation, 
because, first of all, the collision integral is not just a constant and depends on time. Second, 
past (future) infinity in EPP (CPP) corresponds to the UV limit of the physical momentum. 
In this limit g p (rj) harmonics behave as plain waves e ipr> . As the result all the rates inside the 
collision integral can be well approximated by 5-functions ensuring energy conservation [6]. The 
latter fact has a clear physical meaning, because the modes whose wavelength is much smaller 
than the cosmological horizon size, practically do not feel the curvature of the space-time. 
Hence, for them the conditions should be similar to those in flat space. 

Thus, the collision integral becomes negligible as pr] — > 00 both in EPP and CPP. That 
is the reason why in the EPP, independently of the type of the harmonics, one can put the 
moment of turning on selfinteractions, r] , to past infinity. Similarly these observations explain 
the log(l/p?7o) behavior of the leading IR contribution in the CPP, when 77 > I/77. 



B. Kinetic equation in the contracting Poincare patch 

In ([3]) we have found the created particle density n p and anomalous average k p under the 
implicit assumption that throughout the time period [770, 77] both n p and k p have been zero inside 
the collision integral. In fact, we have been using the tree-level propagators in the calculation 
of Q. For the proper consideration of the phenomenon one has to redo the calculation and 
take into account the change of n p and k p in time. If we care only about IR contributions this is 
done with the use of the following ansatz for the IR behavior of the exact Keldysh propagator 

4]: 



D p (r)i,r) 2 ) = (771772) ° 2 1 [hipnx) h*{pr) 2 ) n p (r]) + h(prn) h(pr) 2 ) K p {rj) + c.c] , 77 = y/rjrrj^, (5) 
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where one has to assume the usual separation of scales in the kinetic problem, i.e. that n p 
and hip are slow functions of 77 in comparison with the harmonics h{jprj). Such a separation of 
scales is quite obvious for the massive fields from the principal series, because their harmonics 
do oscillate in the IR limit. However, for the light fields from the complementary series the 
separation of scales is not so transparent, because their harmonics do not oscillate in the same 
limit. 

It is not hard to see that after the assumption about the IR behavior of the exact propagator 
(JSJ), the problem of the calculation of n p and n p reduces to the solution of the Dyson- Schwinger 
equation in the Schwinger-Keldysh diagrammatic technique |4|. In the IR limit, which is of our 
interest, the Dyson-Schwinger equation reduces to a variant of the kinetic equation. One can 
write the kinetic equation both in EPP and CPP % la] . With the use of the following matrixes: 



N p (r] U 772) = % 



D-l 
2 





'0 




)•'-! 




: 









(6) 



n p (r)i)h*(pr} 2 ) «p(ffr) Kp^) 
Kp{vi) h{PV2) M771) h *(PV2, 
the system of kinetic equations for n p and k p can be written compactly. The real equation has 
the form: 



dn p {rj) 



2 A 



d D - l q 



drj' 



Re Tr 



d V J (2vr)(^) J Q (jfrf) 

pM*-M [{ l + N l) N q N p _ q - N* p (1 + N g ) (l + N p . q )] (77,77')+ 
+2 C q ^ p>p ( V ) [N* q (1 + iV ? _ p ) (1 + N p ) - (1 + N* q ) N q . p N p ] ( V , i) + 
+D p , q , p+q (v) [(1 + N P )(1 + N q ) (1 + N p+q ) - N p N q N p+q ] (77, V ') } + [N -> P N] . (7) 
At the same time the complex equation is as follows: 



dK p (rj) 
dr) 



-2 A 



dP~\ 



(2tt)( d - 1 ) J (jtrf) 



a 



(v) 



[l + N p ) {1 + N q ) (l + N p _ q ) 



v drf 
~d KP ->■ -P) Tr 



N N N 



+2 C* qtq _ p>p (v) [N* q (1 + N q _ p ) (l + Np) - (1 + N*) N q . p N p \ (77, 7/) + 

+ D * P , q , P+q (v) [(1 + Np) N* q N* p+q - N p (1 + N*) (1 + N* p+q )] ( V , i) } - [N P N] . (8) 

The term [N — > P N] means that we have to add to the explicitly written expressions in (JTJ) and 
OHJ) the same quantities with every iV substituted by the product P N; Re Tr means that one 
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has to take the real part and the trace of the expression following after these signs; (p — > —p) Tr 
means that one has to take the trace and add to the expression following after these signs the 
same term with the exchange p — >■ —p. Finally Ck 1 k 2 k 3 (v) = ' h*{kiv) hik^f]) h(k 3 r]) and 

Dkrkikaiv) = V^ 21 Hhv) h(k 2 rj) h(k 3 r]). 

Actually the system (J7J) and (jHJ) represents just a preliminary version of the kinetic equa- 
tions. It is not yet suitable to sum only leading IR terms in all loops. It contains the latter 
contributions, but on top of that it also accounts for some of the subleading expressions. In 
particular, one should restrict the dq integrals in (jJJ) and (EJ) to the region q S> p. It is this 
region from where the leading IR contribution does come. For more details see {4, 6] (see also 
the discussion in the next subsection). 

Due to the presence of k p , in ((Tj) we have extra terms on top of those which are shown in 
(j3J). All of them can be obtained from (j4]) via the simultaneous substitutions of (1 + n k )'s and 
n k s by /tfc's or k^s. E.g. we encounter terms of the type: 



[(1 +n p )K k n p ^ k -n p K k (l + n p - k )] . (9) 

The meaning of is that it describes the following two competing processes — the particle 
with the momentum p is lost (gained) in such a situation, in which instead of the creation 
(annihilation) of the two particles, k and p — k, we obtain single p — k excitation and the 
missing momentum k is gone into (taken from) the background quantum state of the theory. 
Similarly we encounter terms of the type 



[(1 + Up) K k K p _ k - n p K k K p _ k ] (10) 

which describe the processes in which both momenta k and p — k are coming from (going to) 
the background state. These observations, in particular, justify the interpretation of k p as the 
measure of the strength of the backreaction on the background state of the theory. 

C. Solving the system of kinetic equations 

Spatially homogeneous kinetic equations in EPP and CPP can be trivially related to each 
other, because to perform the map from EPP to CPP one just has to flip the limits of drj' 
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integration inside the collision integrals in (J7J) and (J5J). On top of that it is necessary to make 
the change of h(x) to h*(x), because the positive energy harmonics are exchanged with the 
negative ones under the flip of the time direction. As the result, solutions in both patches also 
can be mapped to each other. That is physically meaningful because the system of kinetic 
equations d7j) and (JSJ) is valid for the spatially homogeneous states and is imposed on the 
quantities which are attributed to the comoving volume. 

Before continuing with the solution of ([7]) and (jHJ) let us stress that these equations are 
not yet unambiguously defined, because one has to specify harmonics h(x) in their collision 
integrals. This issue is related to the question of the proper definition of quasi-particles in the 
IR limit. One should keep in mind, however, that the kinetic equations and their solutions for 
different choices of h(pq) can be mapped to each other via Bogolubov transformations. 

Intuition from condensed matter should tell one that the proper choice of h(pn) has to 
correspond to the situation in which there is a solution with k p = 0. It was shown in [6( that in 
EPP there is such a solution of the IR limit of © and (IE]), if one chooses out-Jost harmonics, 

n 

h(x) ~ Jin(x). Moreover it is stable under small linearized perturbation of k p [4J. The systems 
of kinetic equations for the other harmonics do not allow such a solution. 

All these observations can be extended to the CPP. Thus, for the out-Jost (in-Jost) harmonics 
in EPP (CPPVthe system (J7J and (jHJ) can be reduced to the single equation ([7]), in which k p is 
set to zero {4] (|: 



x< Re 



IT = fo ND-i / dwiciV) 2 / drjq{qr]) 2 x 

dlogrj (2tt) d V Jo Jo 

{qrjy^Viqrj) (qvT V*(qr)')} \(l + n p )n q - n p (1 + n q ) 2 } (77) + 



+2 Re 



(qrj)^W(qr]) (qr)')~ ifl W(qrf) 



n (1 + n a ) (1 + n v 



{l + n q )n q n p 



(77) + 



+Re 



iffnTViffq) (qvT" V*(q V ') (l + n q ) 2 (l + n p ) - n 2 q n p ( v )}. (11) 



To obtain (TTT1) from (J7j) we have neglected p with respect to q under the integral and Taylor 
expanded h(pr]) and h(pr)') in the limit pr],pr]' — > 0. The latter approximations allow one to 
obtain the equation which accounts for only leading IR contributions in all loops jfj]. Further- 
more, we have subtracted from h 2 (x), (h*(x)) 2 and |/i(a;)| 2 some polynomials of l/\x\ to obtain 



V(x) = h 2 {x) - - ■ ■ ■ and W(x) = \h(x)\' 



const 



.... That makes the collision integral 
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well denned and does not change the IR behavior {(| . 

Eq. ffTTj) does not have the Plankian distribution as its stationary solution, because there 
is no energy conservation. However, one can find a distribution which annihilates the collision 
integral of ( fill) in a different manner: while Plankian distribution in ordinary kinetic equation 
annihilates each term, describing aforementioned couples of competing processes, separately, in 
the present case we will see the cancelation between different terms. The situation is somewhat 
similar to the one in which appears the stationary Kolmogorov scaling solution in the turbulent 
hydrodynamics, as is natural to expect in a system with energy pumping. 

For very small densities the only two terms, which survive in the collision integral of ( TTTT) 
are those, which are responsible for the particle production from the vacuum and for the 

n, 

decay of the produced particles into two [6|. As the result in the EPP the behavior of the 
solution in the vicinity of the stationary distribution (i.e. as prj <C 1) is as follows [6[: n p (ri) = 
C (pi]) ri + 1 , where C is some real constant which depends on the initial conditions; T2 is 
the total particle production rate and Ti is the total particle decay rate. This solution is self 
consistent only for heavy enough particles m ^> 1, because then the equilibrium distribution, 
T2/T1 rj e~ 2nfl , is much less than unity. 

Such a EPP's stationary solution is mapped to 



1 1 



v 



(12) 



after the conversion into the CPP. For the given comoving momentum p it is valid until 77* ~ 1/p. 
As pr] becomes comparable to unity, the approximation, in which ( TTTT) is derived, brakes down. 
But as we have explained above, right after 77* the particle production and decay processes stop, 
i.e. n p remains constant. One encounters a similar situation in EPP: before 77* ~ 1/p we can 
not use fTTTj) . but also before this moment of time the particle production is suppressed. 

In (fl2|) we see the example of the restoration of dS isometry in CPP after the summation of 
all leading loop contributions. The question is under what type of initial conditions the system 
in question drops into such a stage as ( fl2l) in CPP (and/or in EPP). To obtain the answer 
on this question one has to perform the Bogolubov rotation to the out (in) harmonics in EPP 
(CPP) and check if the resulting values of n p and n p are small or not, right at the moment 
when they start to change (which is 770 in CPP and 77* in EPP)^. If they are small enough, then 



t Such a rotation at past infinity of EPP would affect the UV behavior in the theory. But it does not change 



16 



it is natural to expect that k p flows to zero and n p flows to Gibbons-Hawking value in the IR 
limit . 

Note that the number density per physical volume, which scales as l/rj D ~ l , decays in EPP 
and grows in CPP. Hence, it may seem that the backreaction from the produced particles 
eventually will become relevant in CPP. That is indeed true in the case of the presence of 
any slight spatially inhomogeneous perturbation. But in the circumstances of the ideal spatial 
homogeneity that is not true, because the presence of the stationary state means that the 
dS isometry is restored. As the result, in terms of the invariant expressions, such as physical 
volume and physical momentum pq, all invariant quantities remain finite. E.g. the total particle 
density per physical volume is equal to N = rj D ~ l J d D ~ x pnp(j)) and does not explode. 

However, we are going to see now that the kinetic equation in question has yet another 
peculiar solution both in CPP and EPP. In [6] it was shown that in CPP there is a solution 
as follows: 71(77) ~ A-fMogr? ~ r log ' ^ or ver ^ ^ ow momen ta, PV 0. It is valid for i] < 
?7* = e const l x2 ^> l; A is a constant which depends on the initial state, and f oc ^. The 
interpretation of this solution is as follows: Due to the production of particles at low momenta 
(both by the gravitational field and by the decays of particles from higher levels) at high enough 
initial densities, sooner or latter we have the explosion of their density per comoving volume. 

This solution is mapped to 



"<"> ~ (13) 

in the EPP. It is also valid for very low momenta, prj — > 0, for rj > 77*, where i]* = e - const / x2 <^ 1. 

It is worth stressing at this moment that the solution (|T3|) is possible because of two specific 
features of dS space [6|]. The first one is the flatness of the IR spectrum in dS space — weak 
dependence of the harmonics on the momenta in the IR limit. And the second feature is 
that background gravitational field makes possible particle creation and decay processes. In 
particular, one can show that the scattering processes in 4 theory, which are allowed even 
in flat space, do not contribute to the analog of (fl3j) : only the particle production and decay 
processes do contribute. 



the conclusions in the IR limit. After all, note that Bogolubov rotation in the low energy effective BCS theory 
does not affect the UV behavior of QED. 
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Thus, even in the EPP we have the explosive growth of particle density for big enough its 
initial values. The explosion happens at finite proper time and expansion of EPP will not save 
the situation. That will cause a modification of the dS metric. 

Here also it is appropriate to ask the question: under what initial conditions the system 
drops into such an explosive stage both in EPP and CPP? We are not yet in a position to give 
a conclusive answer on this question. 



III. GLOBAL DE SITTER SPACE-TIME 



To begin with we would like to show, following [5] , that in global dS the moment of turning 
on selfinteractions cannot be taken to past infinity. To find the IR divergences at one loop, we 
can use the metric of either fig. 1 or fig. 2. However, to define the kinetic equation in global 
dS we prefer to work with the slicing depicted in fig. 2. Because, as we explain above, for the 
metric of fig. 1 the collision integral vanishes for the times > 1/p and i] + > 1/p, which 
confuses the whole picture of particle kinetics. 

Unlike higher dimensional dS with the spatial slicing of fig. 2, in 2D global dS space the 
one-loop calculation can be made quite explicit. So let us present it before continuing with the 
IR divergences in general dimensions. 

Due to spatial homogeneity of 2D global dS space together with the initial state, we find 
it convenient to perform the Fourier transform of all quantities along the spatial direction 
y? e [0,2tt]: D^ A {t u t 2 ) = f^dife ip ^D K ' R ' A (t 1 ,if;t 2 ,0). Then the one-loop contribution to 
is as follows: 



Di p (ti,t 2 ) = —A 2 / / dt 3 dt 4 cosh(t 3 ) cosh(£ 4 



q= — oo 



Dop{tuh) Do q (t 3 ,t4) D^ p _ q) (t 3 ,tij D 0p (t 4 ,t 2 ) + 



+2 D*(ti, t 3 ) D*(t 3 , U) £$p_,)(* 3 , U) D*(U, t 2 ) + 2 Dg(t u t 3 ) D«(t 3 , f 4 ) D A p _ q) (t 3 , t 4 ) D A (U, t 2 ) + 

I J A A A 

~7 D 0p (tl,t 3 ) D 0q (t 3 ,t 4 ) £> (p-g)(*3,*4) D 0p (t 4 ,t 2 ) - - D 0p (tl,t 3 ) D 0q (t 3 ,U) Do( p - q )(t 3 ,t4) D 0p (t4,t 2 ) (14) 

Here the Keldysh propagator is D^ p (ti,t 2 ) = Re[g p (ti)g*(t 2 )\, retarded and ad- 
vanced propagators are D^ > (ti,t 2 ) = 6 {t x — 1 2 ) 21m [g p (ti)g*(t 2 )] and D p (ti,t 2 ) = 
—0 (t 2 — ti) 21m [g p (ti)g*(t 2 )~\; g p (t) is a solution of the Klein-Gordon equation, which we do 
not specify so far; to is the moment when selfinteractions are turned on. 
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The above one-loop expression can be rewritten as: 



^2 _^ r r+oo 

D? p (h,t 2 ) = -— g p {t x ) g* p (t 2 ) J2jJ t dt * dti cosh (*s) cosh(t 4 )^(t 3 ) 9 P (U) X 



x<^ -e{t x -u)e{t 2 -u) 



9q{h) 9q(U) 9\ P -g\(t3) 9\p- q \(U) + c. 



+ 



o{ti - h) e(t 3 - u) + e(t 2 - u) e(u - h) 



9 q (h) 9* q {U) 9\ P -q\{h) 9\ P - q \(U) - c.c. 



A 2 r r+oo 

ir 9p(h) 9 P (t 2 ) Y] dt 3 dt 4 cosh(t 3 ) cosh(t 4 ) g*(t 3 ) g*(t A ) x 

1 a JJ t0 



x{9(ti-t 3 )9(t 2 -t 4 ) 



9q{h) 9q(U) 9\p-q\(h) 9\p-q\{U) + C 



+ 



9(h - t 3 ) 9{t 3 - U) - 9{t 2 - U) 9{U - h) 



9 q {h) 9* q {h) 9\ P -q\{h) 9\p- q \{U) - c.c 



+ c.c(15) 



We need to find the leading divergent contribution to this quantity in the limit t = tl+t2 — > +00, 
t = ti — 1 2 = const and t — > —00. First, to do that we can substitute tx t2 by t in the arguments 
of the ^-functions. Second, if t — >• +00 and t — > —00, then the leading IR contributions to Df p 
can appear only from the regions of integration over t 3j4 in the vicinity of t and to- 

As we have explained in the introductory section, for arbitrary dimension in these regions 
g p (t) can be approximated by 



9 P (t) 



V h + (pr] + ), 77+ 
r]J 2 h^(pi] ), i]_ 



t ->■ +00 
t —> —00 



(16) 



for the large enough p. Then as pi] — > and pr/ — > the leading IR contribution to Df p can 
be estimated in any dimension. It is contained in the expressions as follows [5|: 
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Dipih, h) = g P (ti) 9p(t 2 ) n p (t) + g^) g p (t 2 ) K p {t) + c.c, 
where n p (t) w ^4§zi ^ rfD_1 <? g^ 4 (r/ 3 7/4)^ x 

xh* + (pvs) K (pr] A ) h* + (qrj 3 ) h + (qr)i) h* + (\q - p\t] 3 ) h+(\p - q\r] A ) + 



' d D 1 q I I dr] 3 r] 4 (7] 3 7] A )^2° x 



D-3 



VO 



(2n) D - 1 

xh*_ (pr) 3 ) (pr] 4 ) h*_(qr) 3 ) h„(qr] A ) h*_(\q - p\r] 3 ) h-(\q - p\r] 4 ), 
and «p(t) « — (27r) £> - 1 7 q J J ^ 4 ^™ 2 X 

x/i + (p»/s) ^+ (to) /i+(gr/s) M? 7 ^) - p\m) h +(\i - p\m) - 

— (27r) £> - 1 J q J V3 J dr]4 ^ 3T] ^ 2 X 
xh*_(pr] 3 ) h*_(pr] 4 ) h*_(qr} 3 ) h_(qr] 4 ) h*_(\q - p\r} 3 ) h-[\q - p\r} 4 ). (17) 

Here 77 = e~* and tjq = e*°. Thus, IR divergence in global dS is just the sum of the contributions 
from the contracting (past infinity) and expanding (future infinity) regions (Poincare patches) 
of the whole space. 

Obviously one can do the same calculation also with the metric of fig. 1 and obtain the 
same as (fTTI) result for the IR divergences, when t — > +00 and to — > —00. 



A. IR divergence for different harmonics 

In this subsection we describe the character of the IR divergences for different h(prj) - 
the solutions of the Bessel equation. That could help to specify the convenient choice of the 
harmonics to specify the proper quasi-particles in the IR limit. 

We start with the Euclidian harmonics which are CPT invariant. I.e. in this case, 

h*_(x) = h+(x) oc Hi^(x), where is the Hankel function of the first kind. Using that 
h + (x) pd A + x lfJl + A_ x~ tfJl as x — >■ (with some known complex constants A±), we obtain the 
IR divergence as follows: 
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n p (t) 



A 2 >Sd-2 



(27T) 



D- 



l log 



p z 7/r/o 



D-3 

dx 3 dxi (x 3 xA 2 x 



j / i 



\A. 



K p (t) 



2 \ 2 Sd-2 , 



o 

X3 
• r 1 



X3 



+ 



•£3 
.r.i 



/i 2 (x 3 ) [^(x 4 )] 2 , 

D-3 

dx^ (x 3 x 4 ) 2 x 
h\x,) [h*(x 3 )} 2 . 



It also follows from the region of integration over dq in (I17p . in which q ^> p. Note that 
to calculate the expression for k v we have used the identity A* + A*_ = A + A-. Thus, for the 
Euclidian vacuum the leading IR divergence is proportional to log (l/r^o) = (t — to) and is 
present both in n p and K p . 

For the in-harmonics h_(x) is proportional to the Bessel function Y ipi (x). Hence, h_(x) ~ 
Ax~ tfl , as x — > (past infinity). At the same time, for this case h + (x) ~ C + x %[1 + C_ x _iAt , as 
x — >■ (future infinity). Here A and C± are some complex constants related to A + and A_ via 
the corresponding Bogolubov rotation. 

Then the leading IR contribution in this case looks as: 



n p (t) 



(2< 



X 



IC 4 



D— 
X3 



I log 

ip 



1 

prj 

\C- 



X3 
X4 



u — o 

dx 3 dx4 (x 3 X4) 2 x 

v + (x 3 )v;(x4) + 



1 ft 



X 2 S 



(2 



7T 



D-2 i 

— log 



,0 



PVo 



K p (t) 



dx 3 dx 4 (X 3 X 4 )^ |A| 2 ( X -± 

\X 3 

2\ 2 S 



ifi 



V*(x 3 )V4x 4 ), 



{2tt) d 

X (x 3 X4 



•>:■>, 



dxd x 



p-a 
2 



x 3 



IfJ 



— I + — 



;r 3 



X4 



v+{x 3 )v;{x 4 ) 



(19) 



where V±(x) = h±(x) 



const 



. {(J. Note that for the in-harmonics n p (t) depends separately 
on t and to rather than on their difference (t — to) and K p (t) does not diverge as to - ► —00. 
Similar situation also appears for the out-harmonics with one crucial difference that now k v 
does diverge as to —> —00, but is finite as t — > +00. One encounters analogous situation for the 
one-loop contribution in the scalar QED on the electric pulse background 15]. 
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For the other a-harmonics the leading IR contributions to n p and k p also, unlike the case of 
the Euclidian harmonics, depend on t and t separately and they both diverge when t — > +00 
and/or when £ — > — 00. 



B. Kinetic equation in global de Sitter 

Having in mind the discussion of CPP in the previous section, it is not hard to derive the 
system of kinetic equations in global dS: 



drij 



{r+co 
Jj,rh-, h ,my,nrh 2 / dt' C^fr) X 



(1 + N*) N n N j2 -N*(l + N n ) (1 + N J2 )\ (t, if) + 
+2 J,,,, lM:jlTl I dt'C* uj2j (t) iV* (l + iV i2 ) (1 + ty) - (1 + iV*) iV i2 N j \(t,t') + 



+00 



dt'D* jJld2 (t) (1 + Nj) (1 + N h ) (1 + N h ) -N 3 N n N n (M') (20) 



and 



dnj 
~~dT 



- [N — )■ PN] — 2 A 2 -™) Tr \ Jj^LT/iifrrto J 

jl,mi;j2,m2 K 



X 



+2 7 - • - 



3 '.myi,"!! ;j2"i2 



d^l jld2 (t)x 

1 + N d ) (1 + N n ) (1 + i\y - iV,- jV^ iV ia ] (t, t') + 
^ Gfi^W [ N h C 1 + (! + ^) - (1 + iVjJ iV,, ty] (t, + 
dt ' D ln,3M + N n N n- N 3 + (1 + N* 2 )](t,t')\(21) 



The terms [N — > PA/] and (m — >■ — m) have been defined in the subsection B of the section II. 
In (T2QT) and (T2~ll) instead of the usual flat space 5-functions, which ensure energy-momentum 
conservation, we have 

^jim 1 \j2m2]j3m 3 = ^ ^^mi (^) ^2*2 (^) X?3™3 (^) > ^3l>32,33 = cos h (^) 9j 1 9j2 9j3i 

1 ' 3\fh\\j2rh2\jzfhz = j ^^jimi (^) ^)2m2 (^) ^j3"i3 (^) ' ^jl,j2,js = cos h (t) (J^ (^j 2 (7j 3 . (22) 
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Unlike the case of EPP and CPP, we do not yet understand what is the appropriate choice of 
the harmonics for this system of kinetic equations. Hence, we are not yet in a position to solve 
it analytically. 



IV. CONCLUSIONS 



We have derived kinetic equations in CPP, EPP and global dS for the massive scalar fields 
from the principal series. Our goal is to see whether the kinetic equation in global dS has any 
stationary solutions. In EPP and CPP we have found stationary solutions. The possibility of 
their existence strongly depends on the ratio between the particle production, r 2 , and particle 
decay, Ti, probabilities. This ratio defines the equilibrium particle density and is proportional 
to the Gibbons- Hawking rate, e~ 2nfl , fj, = \J m 2 — (-^y^) 2 - It is suppressed when the particle 
mass is very big. However, as the mass of the particle m approaches the (D — l)/2 bound, 
the solution is not valid anymore. Furthermore, when such a ratio is grater than one, then the 
solution in question is unstable under the linearized perturbations (4). Thus e.g. for the scalar 
fields from the complementary series the picture of the particle kinetics can be quite different 
from what we observe for the principal series. 

Furthermore, in global dS particle production and decay rates are different from those in 
CPP and EPP. As the result the stationary solutions of EPP and CPP found above are not 
applicable in global dS. The question is if ( 1201) together with ( 12TT) have any other stationary 
solutions. In this respect it is worth recalling that solutions of these equations for the different 
choices of the harmonics are related to each other via Bogolyubov transformations. Hence, if 
the stationary solution does exist, then it means that there is a choice of gj(t) harmonics for 
which the system of equations under consideration has a solution with k p = 0. However, we do 
not see such a choice so far. 

The last, but not least, subject which we would like to mention is the UV/IR mixing. In Q 
it was noticed that UV and IR renormalizations are mixed in the coordinate space correlation 
functions in global dS. In CPP the situation is similar. In fact, the sufficient condition for this 
mixing is the presence of IR divergences of the type discussed above. However, the necessary 
condition for such a mixing is the presence of the non-vanishing collision integral. E.g. one can 
see that in CPP for the case of the stationary solution the mixing disappears and the theory 
becomes renormalizable. 
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We conclude that, in the circumstances when the stationary situation can not be reached 
over the fixed dS background, the resolution of the problem comes only when one takes into 
account backreaction on the dS metric. Because then the stationary stage should be achieved 
in the deformed space. 
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